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We present a novel ground state that explain the continuous modulated charge diagonal order 
recently observed in manganese oxides, at hole densities x larger than one half. In this diagonal 
phase the charge is modulated with a predominant Fourier component inversely proportional to 
1 — X. Magnetically this state consist of antiferromagnetic coupled zig-zag chains. For a wide range 
of relevant physical parameters as electron-phonon coupling, antiferromagnetic interaction between 
Mn ions and on-site Coulomb repulsion, the diagonal phase is the ground state of the system. The 
diagonal phase is favored by the modulation of the hopping amplitude along the zig-zag chains, and 
it is stabilized with respect to the one dimensional straight chain by the electron phonon coupling. 
For realistic estimation of the physical parameters, the diagonal modulation of the electron density 
is only a small fraction of the average charge, a modulation much smaller than the obtained by 
distributing Mn"*"^ and Mn"*"* ions. We discuss also the spin and orbital structure properties of this 
new diagonal phase. 



PACS numbers: 75.47.Gk,75.10.-b. 75.30Kz, 75.50.Ee. 



Oxides of type (Ri_a;Aa;)Mn03 where R denotes rare 
earth ions and A is a divalent alkahne ion, are called 
generically manganites. In these compounds x coin- 
cides with the concentration of holes moving in the Cg 
orbital band of the Mn ions that ideally form a cu- 
bic structure of lattice parameter a. The ground state 
(GS) properties of manganites are determined by the 
competition between at least four independent energy 
scales; the anti-ferromagnetic (AF) interaction between 
the Mn spins, the electron phonon coupling, the elec- 
tronic repulsion and the kinetic energy of the carri- 
ers. In manganites the energy magnitude of these ef- 
fects is the same and very different states can have very 
similar energies. Consequentially by varying slightly 
parameters as carrier concentration, strain, disorder, 
temperature etc, different GS as ferromagnetic metallic 
phases 0, AF Mott insulatorQ, charge and orbital or- 
dered stripe phases 0, 0, 0, |^, or ferromagnetic charge 
ordered phasesQ can be experimentally observed. 

Studies based in pure doubl&^xchange (DE) modellM 
1^, electron-phonon coupling 10], orbital ordering|llj. 
strain interactions 12], or ab initio calculations have 
been successfully applied for explaining the existence of 
ferromagnetic metallic phases at x ^ 0.2-0.4 and the 
existence of the CE phase at half doping. Also for 
large enough AF coupling between the Mn spins or- 
bital orde ring could induce stripes at x—l/m being m 
an integer^J^. In the CE phase the x-y layers are 
AF coupled and into the plains the magnetic structure is 
formed by AF coupled zig-zag chains. The horizontal (x) 
and vertical {y) steps of the zig-zag chains contain three 
Mn ions. In the CE phase the charge is stacked in the 
z-direction whereas in the x-y planes it is ordered in a 
checkerboard form. 

For X >l/2 microscopic models [ill Il6j predict the 
existence of a C phase formed by vertical one dimen- 



sional chains coupled AF. However experiments reveal 
the existence of a diagonal charge modulation, that was 
interpreted!^ 0] as a mixture of two commensurate ad- 
jacent integer period sub-units. On the other hand, by 
similitude with ferroelectric it has been argue that for 
X >l/2 the GS could be formed by commensurate re- 
gions of density x={l — ^) separated by solitons where the 
extra charge is localised 17]. In this picture the stripes 
are formed by an array of solitons. However recently, 
for doping x >0.5, low temperature electron microscopy 
experiments!^ shown an uniform periodicity propor- 
tional to 1 — a;, ruling out the existence of a mixture 
of integer period units joined by discommensurations or 
solitons. Also at x=2/3 a diagonal supercell structure has 
been proposed for modelling neutron scattering data^^ ■ 

In this work we propose a new GS for the mangan- 
ites at doping x—{l — ^) >0.5. In this phase the charge 
is stacked in layers and at each plane it is modulated 
diagonally with only a predominant Fourier component 
(^, ^)^- Magnetically this diagonal phase consists of 
AF coupled x-y layers formed of zigzag chains also cou- 
pled AF. The horizontal and vertical steps of the chain 
are formed by (m + 1) Mn ions. By varying the values of 
n and m this phase describe a GS where the periodicity 
of the charge changes continuously with x. 

Model. In ideal manganites the Mn ions form a cu- 
bic lattice. The crystal field split the Mn d levels into 
an occupied strongly localised triplet and a dou- 
blet of Bg symmetry. Coulomb repulsion prevents dou- 
ble occupancy and aligns the spins of the d orbitals. At 
x=l the Cg orbitals are empty and the superexchange 
coupling between the Mn spins produces an AF GS. 
At a: 7^1 the hopping of electrons between empty eg 
states is possible. The Hund's coupling between the 
spins of the carriers and each core spin is much larger 
than the kinetic energy, and each electron is forced to 
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align with the core spin texture. Then the carriers can 
be treated as spinless particles and the hopping ampli- 
tude between two Mn ions is modulated by the factor 
/i 2 = COS 2^ COS 2^+e^(<^i-<^2) sin ^ gin ^ where 0J 
are the Euler angles of the, assumed classical, Mn core 
spins {Sj} . This is the so called DE model 

For obtaining the GS of the system we study a 
three dimensional DE model coupled to Jahn- Teller (JT) 
phonons. We also include the AF coupling between the 
Mn core spins Jaf and a Hubbard term U for describing 
the strong inter-orbital Coulomb interaction. 

i.j,a,a' i,a,a' 

+ Jaf S.S, + i^(/3Q?, + Q2^ + Qi,) 

<i.j> i 

+ A ^ {QliPi + Q2iTxi + QziTzi) , (1) 
i 

here Cf^ creates an electron in the Mn ions located at 
site i in the eg orbital a (a=l,2 with l=\x'^ — y"^ > 
and 2=|3z^ — >). The hopping amplitude is finite 
for next neighbors Mn and depends both on the type 
of orbital involved and on the direction u between sites 
i and j {tl[l'> = ±V3tl'-f = ±V34':!^ = itl^^!^ = t 
and 2 = 4/3i with t^ ^ = ^i.i = ^1,2 = 0)^10] .In 
the rest of the paper t is taken as the energy unit. 
In the second term of Eq.lQJ ni^a = Cf^Ci^a- The 
fifth term couples the Cg electrons with the three ac- 
tive MnOg octahedra distortions: the breathing mode 
Qii, and the JT modes Q2i and Q^i that have sym- 
metry x^-y^ and iz^-r^ respectively. Qu couples with 
the charge at site i, pi = CfjJi^a whereas Q2i and 
Q^i couple with the x and z orbital pseudospin density , 
Txi = Cf^Ci2 + Cf^Cii and r^^ = Cf^Cn - C%Ci2., respec- 
tively. A is the electron-phonon coupling constant. The 
forth term is the elastic energy of the octahedra distor- 
tions, being /? >2 the spring constant ratio for breathing 
and JT- modes 20J . In the perovskite structures the oxy- 
gens are shared by neighbors MnOg octahedra and the 
Q's distortions are not independent, being cooperative 
effects very important |2lj. In order to consider this col- 
lective effect we consider the position of the oxygen atoms 
as the independent variables of the JT distortions. 

For a particular electron density, a given set of values 
A, Jaf and ?7, and a texture of core spins {S^}, we solve 
self-consistently the mean field version of Hamiltonian 
and obtain the energy, the local charges {pi} and 
orbital pseudospin orientation {xri, Tzi\. We have solved 
the Hamiltonian for x < 0.5 and with the appropriated 
parameters we have recovered the GS obtained previously 
by different groups H [H Q El HI . Therefore we 
are in conditions for studying the different phases that 
can appear for x > 0.5. 

Results. We have studied different magnetic tex- 
ture states, as the three dimensional (3D) ferromag- 
netic (FM) phase, the A phase formed by FM two di- 
mensional (2D) planes coupled AF, the C and the CE 
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FIG. 1: (Color online) Phase diagram for x—2/3 for the three- 
dimensional DE two orbital model with cooperative Jahn 
Teller phonons. The on-site Coulomb term U is take as zero. 
The symbols CO(CD) and OO(OD) stand for charge and or- 
bital ordered(disordered) respectively. 

phases, the island pha ses [ol l2^ of the form (tt/S, tt, tt), 
the CxEi-x phases^, the Skyrmion phasesQ and the 
3D AF phases. Neither of the GS's obtained by compar- 
ing the energy of these magnetic phases present a diago- 
nal modulation of the charge, in any range of the param- 
eters {A, Jaf, U}. Neither of them present a continuous 
variation of the period of the charge modulation when 
the hole concentration varies. 

In this work we propose, for a;=(l — ■^), a new magnetic 
phase consisting of AF ordered layers composed of AF 
coupled zig-zag chains, with vertical and horizontal steps 
formed by m+1 Mn ions. As the tunnelling amplitude de- 
pends on the direction between sites {tf 2— 2) i'^ this 
phase the hopping amplitude is modulated with period 
2m along the zig-zag chains, and this make possible the 
existence of Fourier components of the charge of the form 
( — , —) — . Even in absence of interactions, the modula- 
tion of the hopping along the chain produces the appear- 
ance of gaps in the energy spectra at hole concentrations 
1 — ^ and the system behaves as a band insulator In 
the case of m=2 this phase coincides with the CE phase 
When the Hubbard term or the electron-phonon coupling 
is included, the gap at the hole density 1 — — is enlarged 
and the charge density modulation with period ^(a, a) 
is privileged, being the charge modulated diagonaUy ba- 
sically with a single Fourier component. 

The unit cell of this phase contains 4mx4m Mn ions 
and this impose a computational restriction to the stud- 
ied hole densities. We have calculated, for a;=2/3 and 
x=3/5, the energy of the diagonal and other established 
phases for different values of A, U and Jaf- By compar- 
ing energies we build the corresponding phase diagrams. 
In Fig.l we plot the A — Jaf phase diagram for x = 2/3 
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FIG. 2: (Color online) Schematic view of the diagonal phase 
in the x-y plane. Open and solid symbols indicates up and 
down t2g spins respectively. 



FIG. 3: Phase diagram (A — [/), for the three-dimensional DE 
two orbital model with cooperative Jahn Teller phonons, at 
hole concentrations x=2/3 and x=3/5. 



and U — 0. At small values of Jaf and A the GS is 
the uniform charge 3D FM phase, this state is the stan- 
dard FM metallic phase that appear in the canonical DE 
model. At small values of Jaf and large values of A 
the GS is an insulating 3D FM phase where the electric 
charge and the orbital are ordered. This charge ordered 
FM phase is the 3D analo g to the 2D charge ordered 
phase reported in reference [23. At moderate values of 
Jaf and in a wide range of values of A there is a large 
window in the phase diagram where the GS is the diago- 
nal phase. At small A the 3D-FM state is separated from 
the diagonal phase by the A phase. For large values of 
Jaf island phases of the form (tt/S, tt) appear 20] and for 
much larger values of Jaf (not shown in Fig.l) the pure 
AF phase wins in energy. Topologically equivalent phase 
diagram is obtained for x = 3/5. The important point 
to note is that the window of X-Jaf values, where the 
diagonal phase is stable includes the values for which the 
CE phase is the GS at x=l/2. Therefore we expect the 
diagonal phase to be the GS for x >l/2. 

In Fig. 2 we plot schematically the spin-orbital ordering 
in the unit cell of the diagonal phase at a;=2/3. Open and 
solid symbols represent up and down Mn core spins. The 
lobes indicates that the electron occupy 3a;^ — or 3y^ — 

orbitals whereas small circles represent the occupation 
of the x'^ — orbital. The charge density is constant 
along diagonal directions marked in Fig. 2 and periodic 
along the perpendicular direction. At x=2/3 the electric 
charge on the Mn ions have only two values ni or 712, 
and the charge disproportionation of the electron charge 
depends on A and U. For realist values of X=0.6t and 
U—At the difference of charge is ni-n2=0.15. This value 
is very different from the extreme value of 1 obtained by 
distributing periodically Mn^+ and Mn''^ ions. 



The phase diagram above presented correspond to 
[7=0. It is known '2^ that at a;=l/2 the U interaction 
destabilizes the CE phase against the C phase. For ana- 
lyzing the estability of the diagonal phase we plot in Fig. 3 
the X-U phase diagram for x=2/3 and x—3/5. In the C 
phase only an orbital — is occupied and the energy 
of this phase is independent of U. Therefore the on site 
Coulomb interaction favors the C phase against the di- 
agonal phase and it is necessary a finite electron phonon 
coupling to stabilize the diagonal phase. However the 
require values of A are moderate and the diagonal phase 
has lower energy than the C phase for realistic estimation 
of this parameter po||. 

The results shown until here indicate the existence of 
a diagonal phase at x=2/3 and x=3/5. As commented 
above the Hubbard term and the electron phonon in- 
teraction favor the existence of a predominant Fourier 
component, ^(^, ^) of the charge density. In Fig. 4(a) 
we plot the Fourier transform of the electrical charge for 
x=l/2, x=2/3 and x=3/5 and X=0.6t and U=U. Only 
the Fourier components ±(^, have values consid- 

erably different from zero. 

There is experimental evidence 03| that the period of 
the charge modulation varies continuously with the hole 
density. The study of a system with a hole concentra- 
tion close to half filling, as x=0.52 or x=0.58, requires 
the use of very large unit cell and it is computationally 
unaccessible. However for zero electron phonon coupling 
the properties of the diagonal phase can be obtained by 
solving one dimensional system with 2m Mn ions per unit 
cell. In this form we study systems with x close to 1/2. 
In Fig. 4(b) we plot the Fourier transform of the elec- 
tron charge for f7=4t, A=0 and hole densities a;=l/2. 
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FIG. 4: Diagonal, (g, q) Fourier component of the electric 
charge in the diagonal phase, (a) corresponds to X—0.6t and 
U^4:t and (b) to A=0 and f/=4t. 



a;=13/25 and a;=29/50. It is noticeable that even in 
absence of electron phonon coupling the charge density 
is diagonally periodic with a prevailing Fourier compo- 
nent ^(^,^)- We expect that the inclusion of elec- 
tron phonon coupling will make more predominant this 
Fourier component. These results indicate that the di- 
agonal phase describes an state where the electric charge 
is ordered diagonally and the periodicity of the charge 
modulation changes continuously as a/(l — x). 

Summarizing, we propose a diagonal phase for man- 
ganites at doping x=(l — In this phase the charge 
is modulated diagonally with a prevalent Fourier com- 
ponent ^(^j^)- Magnetically this phase consists of 
zig-zag chains with horizontal and vertical steps formed 
by m+1 Mn ions. We obtain that this phase has lower 
energy than the previous suggested phases and is the 
ground state for a wide range of the relevant physical 
parameters. In agreement with recent experimental re- 
sults this phase describes a state where the periodicity of 
the charge changes continuously with the hole doping. 
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